Comparison of s- and d-wave gap symmetry on nonequilibrium superconductivity 
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Recent application of ultrafast pump/probe optical techniques to superconductors has renewed 
interest in nonequilibrium superconductivity and the predictions that would be available for novel 
superconductors, such as the high T c cuprates. We have re-examined two of the classical models 
which have been used in the past to interpret nonequilibrium experiments with some success: the fj,* 
model of Owen and Scalapino and the T* model of Parker. Predictions depend on pairing symmetry. 
For instance, the gap suppression due to the excess quasiparticle density n in the fi* model, varies 
as n 3 / 2 in d-wave as opposed to n for s-wave. Finally, we consider these models in the context 
of SIN tunneling and optical excitation experiments. While we confirm that recent pump/probe 
experiments in YBCO, as presently interpreted, are in conflict with d-wave pairing, we refute the 
further claim that they agree with s-wave. 

PACS numbers: 74.40.+k,74.72.-h,74.25.Gz 



I. INTRODUCTION 

The field of nonequilibrium superconductivity was very 
active throughout the late 1970's to mid 1980's when 
it was realized that novel effects in the superconducting 
state could be induced by converting the electron dis- 
tribution function into a nonequilibrium one..lj Differ- 
ent experimental techniques were used to prepare such a 
nonequilibrium state, for example, tunnel injection and 
optical irradiation, and a body of work arose from both 
experimental and theoretical efforts in this area. A use- 
ful summary of this work near the end of this period of 
time can be found in a book edited by Langenberg and 
Larkinp] and other broad-based texts have also appeared 
more recently |2|, |3j- 

The advent of high T c cuprate superconductivity in 
1986 interrupted work in this and other areas as the 
community turned its attention to this new challenge 
and, consequently, extensive work in the area of nonequi- 
librium superconductivity has languished until more re- 
cently. However, during the period following the orig- 
inal burst of activity, new state-of-the-art experimen- 
tal probes have been developed which provide excellent 
opportunities for renewed interest in this field, not to 
mention the potential for new insights provided by the 
new generation of materials exhibiting novel supercon- 
ductivity, such as the cuprates. Some of these probes 
which can be turned to this problem are: STM, ultra- 
fast lasers, spin-polarized tunneling injection, terahertz 
spectroscopy, etc. 

As early as the mid-eighties, the pump/probe fem- 
tosecond spectroscopy was exhibiting its potential as a 
technique for investigating nonequilibrium phenomena in 
metals and superconductors. In these experiments, an ul- 
trafast laser pulse (~ 100 fs) incident on a sample as a 
high energy "pump" quickly excites the electrons out of 
equilibrium which then relax back to thermal equilibrium 



with the lattice via the electron-phonon interaction. An- 
other laser pulse delayed in time "probes" the system of 
electrons by reflection or transmission spectroscopy. As 
the system of electrons relaxes, the transient reflectivity 
or transmissivity decays with time over a scale of picosec- 
onds or less allowing this experiment to probe carrier dy- 
namics in a time-resolved fashion. A theory was proposed 
by P.B. Allen|j] for the relaxation of quasiparticles in the 
normal state, which could be measured in these experi- 
ments, resulting in the extraction of the electron-phonon 
renormalization parameter A (as the quasiparticles relax 
through interactions with the system of phonons). Ex- 
periments were performed which measured this param- 
eter using Allen's theory and excellent agreement was 
found with other values in the literature for both ordinary 
metals and superconductors in the normal state. Q In- 
deed, this parameter was measured for the first time in Cr 
by this technique. This extraordinary success has led 
experimentalists to use the femtosecond laser as a probe 
of high temperature superconductivity |6l 171 Id lijl. I U| and 
in general several groups have been developing ultrafast 
techniques of similar sort for measuring nonequilibrium 
phenomena in superconductors |llUl2fl . 

Here, we are interested in the state that arises when 
the nonequilibrium excitations, created by a laser pulse 
or by tunneling injection, have fallen to the gap edge but 
have not yet recombined into the condensate (bottleneck 
effect). In the first case, there is some debate amongst 
experimentalists as to whether the high energy laser used 
for pumping and probing can truly measure the distribu- 
tion of quasiparticles at low energy and several groups are 
developing techniques to probe at lower energy of order 
of the gap to address this issue. 

The main thrust of our work has involved the use of two 
models employed in the past to describe a nonequilibrium 
distribution of quasiparticles: the T* model of Parker |13j 
which uses an equilibrium distribution function at an ef- 
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fective temperature T* relative to the bath temperature 
T and the tz* model, originally proposed by Owen and 
ScalapinojlJ], where the system is described in terms of 
a new chemical potential for the excited quasiparticlcs. 
The former approach has been used by Kabanov et a/.|3 
to analyze their optical data, whereas the latter approach 
has been used for systems where excess particles are in- 
jected into tunnel j unctions [l5j. While these two models 
are somewhat simplified, they appear to have been effec- 
tive in capturing some of the experimental results on low 
T c superconductors. 

In section II, we calculate how the superconductivity is 
modified as a function of the nonequilibrium excess quasi- 
particle number density n. This leads to modifications in 
the gap which we calculate numerically for various values 
of temperature T characterizing the sample before irra- 
diation as a function of n in both ji* and T* models and 
for s- and d-wave. For T — and in the limit of n — > 0, 
we also obtain analytic results for the gap reduction ver- 
sus n, for the chemical potential in the \i* model and 
for the nonequilibrium effective temperature for the T* 
model, as well as for the free energy difference between 
the nonequilibrium superconducting state and the corre- 
sponding equilibrium normal state. The analytic limits 
are tested against the numerical work and found to be 
close to the exact results even as n increases towards 
its critical value where superconductivity is destroyed. 
Results for d-wave are compared with s-wave and im- 
portant differences are established. In section III, as an 
explicit example of an application of our results, we con- 
sider a S-I-N tunneling junction with a nonequilibrium 
state on the superconducting side which we assume can 
be described by a fi* model. We show that the current 
voltage characteristics are modified in two major ways. 
First the amplitude of the gap is reduced because of the 
presence of a nonequilibrium number of excess quasipar- 
ticles n and secondly the entire characteristic is shifted 
upward by a factor of n in appropriate units. Also the 
voltage at which the current is zero can be used to mea- 
sure the chemical potential /i*. Separate measurements 
of the gap reduction, the chemical potential, and the up- 
ward shift in I-V characteristic would allow a consistency 
test of the model. In section IV, we consider the specific 
case of pump/probe experiments and agree with previous 
theoretical workjjj that the existing data, as currently in- 
terpreted, is not consistent with d-wave gap symmetry, 
but disagree that it is consistent with s-wave. In a final 
section V, we draw conclusions and give a summary of 
our results. 



II. THEORY 

We consider two models used in the past for the treat- 
ment of non-equilibrium superconductivity. For an s- 
wave BCS superconductor, Owen and Scalapino consid- 
ered a state in which there exists a finite distribution of 
excess quasiparticles at the gap energy in addition to a 



condensate. In their //* model, [lj] thermal equilibrium 
is assumed although chemical equilibrium is not for the 
paired and unpaired electrons. This is mimicked through 
the introduction of a chemical potential /x* in the Fermi 
function which represents a constraint on the quasiparti- 
cle excitation number. With this chemical potential the 
Fermi function is: 

f(E k - f i*) = [l + cxpP(E k -^)}- 1 (1) 

with the BCS gap equation modified to be: 

' r tanh(/3(S fe - /i*)/2) (2) 



N(0)V J y/ei + A^n) 

where V is the pairing potential, N(0) is the electronic 
density of states at the Fermi surface in the normal state 
and the excess quasiparticle density n is given as: 
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A(0) Jo 



[f(E k - /i*) - f(E k )]de k 



(3) 



where f3 = l/(k B T), E k = y 1 e\ + A 2 (n), and k B is the 
Boltzmann constant. Here n is measured in units of 
4A(0)A(0). The 4 is introduced for spin and for particle- 
hole parts of the excitation spectrum. A(0) = A(n = 0) 
is the superconducting gap in the equilibrium state, fi- 
nite and isotropic over the entire Fermi surface for s-wave 
gap symmetry. This model will be applied later to discuss 
tunneling. 

Alternatively, Parker^^ considered a T* model where 
instead of a /i* in the Fermi function, a T* is used: 



f(E k ,T*) = [l+eME k /k B T*)}- 



(4) 



with the other equations modified accordingly. This 
model is the one used by Kabanov et ai.|2| in their anal- 
ysis of the pump/probe data. 

We consider first, the n* model for an s-wave BCS 
superconductor. At zero temperature the existence of the 
excess quasiparticles perturb the condensate by blocking 
states which would otherwise be available to form the 
condensate in a variational sense, and this lowers the 
value of the gap. The exact gap equation and relationship 
between chemical potential and n are respectively, 



A(n) _ 

A(0p I A(0) 



1 and nA(0) = \J ii* 2 - A 2 (n) 

The first equation in (5) comes directly from the gap 
equation (2) evaluated at zero temperature with reference 
made to the equilibrium case which allows us to elimi- 
nate the pairing potential in favour of A(0). The second 
follows from Eqn. (3). The grand potential s (n) (the 
familiar formula is given later for the anisotropic case in 
Eqn. (13)) in the isotropic case (at T = 0) is 

"TTTTvT = ^777v^ = {n)-2fi vV ~ A 2 (n), 



iV(O) 



iV(O) 



(6) 
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TABLE I: Analytical forms for n -> at T = in the Li* 
model. Note n is in units of 4AT(0)A(0), where N(0) is the 
single spin density states and A(0) is the T = and n — 
gap (maximum in d-wave). 

fx* model s-wave d-wave 



A(n)/A(0) 
2AF(n)/iV(0)A 2 (0) 
M */A(0) 



s-wave 
l-2n 
-1 +8n 
l-2n 



4^3/2 
1 + ^ n 3/2 



1 - 
+ 



^2: 



1/2 



where this is the difference between the nonequilibrium 
superconducting state and its normal equilibrium coun- 
terpart (i.e. with no excess quasiparticles). The differ- 
ence normalized to the equilibrium superconducting state 
condensation energy is 



2AO(n) _ 2[n s (n) - fl N (n = 0)] 



Af(0)A 2 (0) 



7V(0)A 2 (0) 



-1 + 8n 2 (7) 



to lowest order in n. To obtain Eqn. (7) we have used 
expressions for A(n)/A(0) and for /x*/A(0) valid to sec- 
ond order in n. They are A(n)/A(0) = 1 — 2n — 2n 2 and 
H*/A(0) = l-2n-3n 2 /2 (entered in Table I to lowest or- 
der). If we add to the grand potential, Af2(n), the num- 
ber of excess quasiparticles multiplied by the chemical 
potential, i.e. /i*n where n is the first term of Eqn. (3), 
normalized in the same way as Eqn. (7), we get the nor- 
malized free energy difference at zero temperature which 
we denote by 2AF(n)/N(0)A 2 (0). This is evaluated to 
be 



2AF(n) 
7V(0)A 2 (0) 



-1 + Sn 



(8) 



(entered in Table I). 

In the top frame of Fig. 1, we present our numerical 
results for the ratio A(n)/A(0) as a function of excess 
quasiparticles n (solid curve) and compare with the ap- 
proximate result A (n)/ A (0) = 1 — 2n (dashed curve). We 
see excellent agreement at small n. As n is increased, the 
continuation of the solid curve is denoted by the dots. It 
is terminated at the point where the free energy for the 
nonequilibrium state becomes equal to its normal state 
value and a first order transition occurs. This can be 
seen more clearly in the bottom frame which shows the 
normalized free energy difference of the nonequilibrium 
(n ^ 0) state, 2AF(n)/N(0)A 2 (0) as a function of n. 
The solid curve applies to the exact result at T = while 
the dashed is the approximate result (Eqn. (7). which fits 
well the exact result at small n and is semiquantitative 
in the entire physical region. The first order phase tran- 
sition to the normal state occurs at n c ~ 0.15. The con- 
tinuation of the solid line for the free energy difference to 
values of excess quasiparticles n beyond the critical value 
is indicated by a dotted curve just as in the top frame 
for the gap. We note that both the gap A(n) and the 
free energy difference AF(n) as a function of n fold back 
on themselves beyond a certain value of n, but that the 
free energy remains positive for the entire dotted region 
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FIG. 1: Top frame: A(ra)/A(0) versus n, at T = for the 
fi* model with an s-wave gap. The solid curve is physical, 
the dotted curve is not. This latter curve represents the case 
where the free energy of the normal state is lower than that 
of the superconducting state as shown in the bottom frame. 
The presence of excess quasiparticles suppresses the gap and 
eventually leads to a first order transition to the normal state 
at n — 0.15. Bottom frame: AF = Fm — Fs, the free energy 
difference, versus n. In both frames, the dashed curve is the 
small n limit (see Table I). 



i.e. the nonequilibrium state has higher free energy than 
does the normal state (AF(n) > 0) in this region. 

Now we treat the d-wave case. The equation relat- 
ing fi* to n is nA(0) = f^* N(E)dE where for small 
E, N(E) ~ E/A(n). Here, A(n) is the maximum d- 
wave gap where the gap A(<j>) at any point <j> (the polar 
angle for momentum) on the 2-dimensional Fermi cir- 
cle in the CuC>2 Brillouin zone is A((f>) = A(n) cos(2<^>)) 
with zeros in the (7r,7r) direction and other symmetry re- 
lated points. The small n limit gives ^*/A(0) = \f2n x ^ 
which differs radically from the s-wave case and reflects 
the gap symmetry with nodes (see Table I). Numerical 
results for /x*/A(0) versus n are given in Fig. 2. The top 
frame applies to the s-wave case and is for comparison 
with the bottom frame for d-wave. The dashed curves in 
both frames are our approximate analytical results which 
are seen to match well the exact results (solid curve for 
T = 0) in the small n limit. The remaining curves are at 
finite temperature T as indicated in the caption, namely 
T/T c = 0.3, 0.5, 0.7 and 0.9. Several features are worth 
noting. For s-wave, the zero temperature behaviour of 
the chemical potential as a function of n is qualitatively 
different from the case for finite temperature. In the 



4 



o 

<T 



71 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r— 
/z'- mo del 




0.8 



0.6 - 



0.4 



0.2 



0.0 



I I I I I I I I I 



0.00 




d-wave 



0.20 



FIG. 2: The parameter fj,* versus n for several temperatures 
shown for the s-wave (top frame) and d-wave (bottom frame) 
gaps. The dashed curve is the small n limit (see Table I). 
From top to bottom, the solid curves are for T/T c = t = 0, 
0.3, 0.5, 0.7, 0.9. Here only the physical part of the curves 
are shown. 



limit of n — > i.e. very few excess quasiparticles, the 
chemical potential must clearly be equal to A(n = 0) at 
T = 0. In this case the lowest energy available quasipar- 
ticle states are at A(0) where there is an inverse square 
root singularity in the density of states and hence all the 
excess quasiparticles can be accommodated at the gap 
energy. As n increases out of zero, the gap A(n) in the 
nonequilibrium state decreases from its value at n = 0. 
The inverse square root singularity shifts to lower en- 
ergy and there are now many states at and around A(n) 
and it turns out that all the excess quasiparticles can be 
accommodated in a small energy range around the new 
gap value. We have already noted that to second or- 
der in n, A(n)/A(0) and /x*/A(0) differ by a factor of 
n 2 /2, specifically, jti*/A(0) = [A(n)/A(0)] + n 2 /2 which 
implies that \i* falls a few percent above the nonequilib- 
rium value of the gap in units of A(0). Note that the 
inequality A(n)/A(0) < /U*/A(0) (at T = only), found 
to hold to second order in n, was also verified in the nu- 
merical work, which shows that the difference between 
A(n) and /1* are always small even outside the validity 
of our expansion. That this difference should be small is 
a reflection of the square root singularity in the density 
of states. 

The situation is very different in the d-wave case and 
in s-wave at finite temperature. In these two cases the 
chemical potential becomes small as n —* 0. For the 



d-wave case this is easily understood because there is a 
small but finite density of states at any nonzero value of 
energy w^O. The excess quasiparticles can occupy these 
states and hence fi* — > as n — > 0. For the s-wave case 
at finite T a different argument holds. In this case the 
thermal factor f(E? — fj,*) gives the probability that the 
state Ej: is occupied at finite T. This probability can be 
increased over its value for ^* = simply by having fi* 
take on a small finite value to accommodate the excess 
quasiparticles. At low temperature, however, the ther- 
mal tails of the occupation factor are small in the region 
of the gap and fx* must increase fairly rapidly as n in- 
creases. This is seen most clearly in the second highest 
curve in the top frame of Fig. 2 which corresponds to 
T/T c = t = 0.3. Also as the temperature is increased 
decreases as expected. In the d-wave case shown in 
the lower frame of Fig. 2, fj,* starts from zero at n = 
even at zero temperature because, as we have already in- 
dicated, there are states available at any energy above 
u> = 0. Comparing top and bottom frame we note that 
the chemical potential for t = 0.3 (to be specific) rises 
more rapidly in the s-wave case and becomes bigger than 
for d-wave. This can be traced to the fact that for d- 
wave the part of the density of states that is occupied 
by the excess quasiparticles is in the range to fi* while 
in the s-wave case it is the region just about the gap 
A(n) which is relevant. As the temperature is increased 
towards T c , the differences in the quasiparticle density 
of states between s- and d-wave become smaller and the 
chemical potentials start to become very similar. A sec- 
ond feature to be noticed is that at finite T the curves for 
fi* extend to higher values of n for the s-wave case than 
they do in the d-wave case although the reverse is true at 
zero temperature. In all cases the curves terminate when 
the free energy difference between normal and nonequi- 
librium superconducting state becomes zero or there are 
two solutions and the one with the lowest free energy is 
chosen. This occurs at smaller values of n for the d-wave 
case as compared with s-wave for the given temperature 
T^0 shown. We will return to this issue later on in our 
discussion of Fig. 4. 

The gap equation with a pairing potential of the form 
V^ck' = V cos(2(j)') cos(20), where k is momentum on the 
Fermi surface, with a distribution of excess quasiparticles 
included through the introduction of a chemical potential 
takes the form 



1 



N(0)V 



/o y/ei + tf{n) cos 2 (20) V * » f 

(9) 

with E k = y/el + A 2 (n) cos 2 (20). The bracket (• • • ) in- 
dicates the angular average and e k is energy integrated 
in a rim of width uj c about the Fermi energy. With ref- 
erence to the n = case (i.e., fj,* — 0) we can rewrite 
Eqn. (2) to read at T = 



VA(0)J 



-4cos 2 (2(/))de 



v/e 2 + A 2 (n) cos 2 (2(f>)/ E k <^» 



(10) 



■5 



where the integration over energy and <f> must duly take 
account of the restriction E < \i* . For small n — > the 
leading order gives (A(n) = A(0) + SA(n)) 



5A(n) 



8 r !A 

^" »/ cos -1 ("m* 

f 

J|A(n 



cos 2 

'cos-^^i'/Afn)) 

dE 



' | A(n) cos 0'| 



^/E 2 - A 2 (n) cos 2 , 



(11) 



where we have changed from <j> to $ — 2<j>. But the 
lower limit in the </>' integration in Eqn. (11) restricts 
the integration to the nodal region which corresponds to 
<t>' = tt/2. We find 



5A{n) 



--( 

7T V 



4V2 



A(n) 



x 2 (ix In 



i + VT 



-n 



3/2 



(12) 



(entered in Table I) where we have used the relationship 
fj,*/A(n) = \[2n}-l' 2 to lowest order. In Fig. 3 we show ex- 
act numerical results for the normalized gap A(n)/A(0) 
as a function of n for the d-wave case (solid curve) 
and compare with our approximate result (dashed curve) 
which applies only at small n. The agreement is excellent 
even up to the point where the first order transition to 
the normal state occurs. This is where the solid curve is 
extended into the dotted curve. The gap function as a 
function of n is reduced less in d-wave (Fig. 3) as com- 
pared to s-wave (Fig. 1) all the way to n = n c . The free 
energy difference AF(n) becomes zero at n = n c ~ 0.17 
which is to be compared with ~ 0.15 in the s-wave case. 
At the critical n, A(n)/A(0) is almost 0.6 for s-wave 
while in the d-wave case it has not yet reached 0.8. The 
blocking of states by the excess quasiparticles has much 
less effect on the condensate wavefunction as reflected in 
the change in the value of the gap in d-wave than in s- 
wave because now the excess quasiparticles accumulate 
in the nodal region. Since the gap is zero or near zero 
in that region, it is clear that these states do not con- 
tribute much to the lowering of energy brought about by 
the formation of Cooper pairs. 

To establish where this first order transition occurs, we 
need the free energy. The formula for the grand potential 
for the superconducting state with n excess quasiparticles 
is 

n s (n) = 2k B TY,Ml-.f(E k -^)) 

k 

+ -E k + £jj-(l- 2f(E k - /i*))](13) 

k 

and for the normal state with n = it is 
n N (0) = 2fc B T^ln(l - /(M)) + 5> fc - \e k \) (14) 
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FIG. 3: The fx* model at T — for a d-wave gap with the 
curves labelled in the same manner as for Figure 1. The 
gap is suppressed less rapidly in d-wave. The presence of 
excess quasiparticles, which normally weaken the condensate 
by blocking states, are less effective in interfering with the 
formation of the superconducting wavefunction in d-wave as 
they accumulate at the nodes, in the first instance, which is 
a region where the gap is close to zero. 



N(Q) taken out of the integration. In the limit n — > 

Ag(n) = n s (n)-n N (0) _1 2 
N(0) ~ N(0) 4 (n> 

+ 4 r N{E){E- n*)dE- i/A 2 (n) (15) 
Jo 2 

where / is the same integral as appears on the right- 
hand side of Eqn. (11). The first term in (15) is the 
usual expression for the condensation energy of a d-wave 
superconductor but with A(n) = A(0)(1 - 4\/2n 3 / 2 /3) 
replacing the gap amplitude A(0) which applies to n = 0. 
In An(n = 0)/N(0) only A 2 (0)/4 enters. The two extra 
terms in Eqn. (15) can be worked out analytically as 
n — > and lead to 



Af2(n) 



.*3 



1 A 2 / \ 2 fj,- 

-4 A (n) -3AR 



If- 

3 VA(n) 



A 2 (n) (16) 



only in the first term on the right-hand side of the 
equation must we retain the n dependence in A(n). 
The difference in grand potential AQ.(n) normalized to 
A 2 (0)iV(0)/4 is easily worked out to be 



The sum over k can be converted to energy and the con- 
stant two dimensional electron density of states factor 



6 



The normalized free energy AF is obtained by adding 
H*n to Eqn. (13) and after normalization we get 



4AF(n) 
iV(0)A 2 (0) 



32^2 



,3/2 



(18) 



which is entered in Table I. Numerical results at any 
value of n are shown in the bottom frame of Fig. 3. The 
solid line is our numerical result for 4AF(n)/7V(0)A 2 (0) 
at T = and the dashed curve our approximate result 
(Eqn. (18)). The analytic result agrees well with the full 
numerical solution at small n and differs slightly near the 
critical value of n = n c where the first order transition 
to the normal state occurs at n r as 0.17. 
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FIG. 4: The ratio of A(n, T)/A(0, 0) versus n for finite tem- 
perature in the n* model. The top frame is for the case of 
an s-wave gap and the bottom frame is for d-wave. Curves 
are shown for T/T c = t = (solid curve), 0.3 (dotted), 0.5 
(short-dashed), 0.7 (long-dashed), 0.9 (dot-dashed). Only the 
physical part of the curves are shown. 

In Fig. 4 we show our numerical results for the gap as 
a function of n at various temperatures. Both frames are 
for the ll* model. The temperatures are T/T c — t = 
(solid curve), 0.3 (dotted), 0.5 (short-dashed), 0.7 (long- 
dashed) and 0.9 (dot-dashed). The top frame is for s- 
wave and is for comparison with the bottom frame which 
is new and applies to d-wave. Note that for s-wave, the 
T = curve is below the dotted curve for t — 0.3. This 
agrees with findings of Owen and Scalapino and has its 
origin in the blocking process referred to previously. At 
zero temperature the excess quasiparticles block impor- 
tant states which cannot be used in the coherent superpo- 
sition of states which form the Cooper pair condensate. 



At finite temperature the blocking is less effective be- 
cause it is the states closest to zero energy that are the 
most effective in forming the condensed pairs while the 
thermal factor depopulates these states. By contrast, for 
d-wave, the T = curve is above the t = 0.3 (dotted 
curve) as we have already noted. In this instance the 
blocking at T = is much less effective and consequently 
temperature is not as important an effect. We note again 
that, at T = 0, the d-wave gap is reduced less than in 
s-wave for the same value of n and that the critical value 
of n, at which a first order transition to the normal state 
takes place, is larger. At the higher temperatures shown, 
however, the reverse holds. Also, note that as the tem- 
perature rises towards T c the difference between s- and 
d-wave get less pronounced as the differences between the 
two quasiparticle density of states become small and also 
more states are involved. 

The nonthermal quasiparticle distribution used in the 
yit*-model has an interesting aspect in that it allows for 
the system to become unstable to quasiparticle density 
fluctuations. [l6L \\% Essentially, if the quasiparticles are 
injected uniformly in the sample, the density fluctua- 
tions will act to draw off quasiparticles from some re- 
gions thereby increasing the superconducting gap locally 
and flowing those quasiparticles to other regions, caus- 
ing an accumulation which lowers the local gap, possibly 
even driving the local region normal. This phase sepa- 
ration could be either a static or a temporal structure. 
Such a state has been studied initially by Chang and 
Scalapino |l6j| and Scalapino and Huberman|l7l| for the s- 
wave superconductor and experimental verification of a 
density instability leading to an inhomogeneous multigap 
state has been done by several groups 18] using tunnel in- 
jection in thin film nonequilibrium superconductors. The 
theoretical signature of such an inhomogeneous state in 
the ^*-model is that dLi*/dn\ T < 0.|H|lZl From Fi g- 2 > 
we find that the variation of [i* with n differs in s- and 
d-wave and by examining the slopes of these curves, in 
particular, the point where the slope goes negative, we 
can reproduce the s-wave phase diagram of Chang and 
Scalapino|l6j|. shown in the upper frame of Fig. 5, and 
provide the equivalent prediction for d-wave in the bot- 
tom frame. The dashed curve in these phase diagrams 
marks the boundary between the normal state (NS) and 
the superconducting state (either homogeneous or inho- 
mogeneous). This boundary is entirely a first-order tran- 
sition. The area labelled IN, is the region of n and T, 
where the slope of the chemical potential curve is neg- 
ative and an inhomogeneous state is predicted to exist. 
The solid line marks the boundary between it and the ho- 
mogeneous superconducting state (SC). There are quali- 
tative differences between the s- and d-wave cases. The 
region of the inhomogeneous phase is quite large in the 
s-wave case and almost non-existent in d-wave and at 
low temperature the s-wave superconductor would likely 
be phase separated whereas, the d-wave one would not 
be. While the inhomogeneous state may be of interest to 
study in itself, in the d-wave case it may be encouraging 
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to note that attempts at experimental verification of our 
predictions for power law dependences, and other results 
presented in this paper, are unlikely to be hampered by 
the presence of an inhomogeneous phase. 
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FIG. 5: The phase diagrams calculated in the the fi* model 
for the s-wave (top) and d-wave (bottom) gaps. Based on the 
slope of fi* versus n one can determine the region of the phase 
diagram where there is a homogeneous (SC) and an inhomo- 
geneous superconducting state (IN). The transition from the 
superconducting state to the normal state (NS) is always first 
order and is represented by the dashed line. 

Next we consider briefly the case of the T* model which 
is just a simple heating model if only the electronic sys- 
tem is considered. Similar approximate analytic calcula- 
tions can be done to get various relationships in the limit 
n — > for the case when the sample before irradiation 
is assumed to be zero. These analytic derivations are 
supplemented with full numerical work in which we also 
consider the case when the sample is initially at finite 
temperature T. 

We begin with the s-wave case and return to the gap 
equation shown in Eqn. (2), now modified according to 
Eqn. (4) rather than Eqn. (1). In the limit of T — > 0, the 
result for the lowest order correction to the gap is well 
known 19] : 



SA{n) 



\ 2-Kk B T* Am/ksT * 



A(0) 



(19) 



The relation between n and T* can be trivially obtained 
as n = ^7rT*/2A(0)e- A (°)/ feBT * and so <5A(ri)/A(0) = 
-In, 



TABLE II: Analytical forms for n -> at T = in the T* 
model. Note n is in units of 4JV(0)A(0), where N(0) is the 
single spin density states and A(0) is the T = and n — 
gap (maximum in d-wave). 
T* model s-wave d-wave 

A(n)/A(0) T - ^ 1- f!(3n) 3/2 

T* vs. n n = O.M^/T*/T c e- 1 - 76T - /T " T* /T*= 2.36n 1/2 



The d-wave case is not as well known and we include 
the critical steps here 



lQl .A(0)J 



tt/2 



2d<f>' 2 
COS 



dEe- E / k * T ' 



/A(ri)cos0' \/E 2 - A 2 (n) COS 2 (f)' 

which can be manipulated into 
'A(n)\ « z-^ 2 



(20) 



In 



A(0)J 



COS 2 < V e -A(n)ea.*7*sT' 



dxe -x/k B T' 



lo y/ x(x + 2 A(n) cos 4>') 
which can be manipulated into 
'A(n)\ « W 2 



(21) 



hi 



A(0) 



COS 2 ^ e ~A(n)cos0'/fc B T* 



die-* 1 " 



o y x(x + 2 A(n) cos 4>') 



(22) 



The integral over <fi' is peaked around cos0' = 0, i.e. 
near 7r/2 which allows us to approximate it by 



In 



A(n) 
A(0) 



dyy 2 e- A ^y/ kBT ' 

dxe~ x / kBT * 
o y/xjx + 2A(n)y) 



from which we get 

6A(n) 

Also from the definition of n we get immediately 

7T 2 / T* \ 2 



(23) 



(24) 



(25) 



Exact numerical results agree well with these approxi- 
mate n — ► expressions which we summarize in Table II. 

In Fig. 6 we show numerical results for A(n, T)/A(0, 0) 
versus n where we have normalized the maximum gap 
A(n, T) to the zero temperature equilibrium case. The 
top frame is for s-wave while the bottom is d-wave. In 
each frame the short dashed curve is the approximate 
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FIG. 6: The ratio of A(n,T)/A(0, 0) versus n for finite tem- 
perature in the T* model. The top frame is for an s-wave 
gap and the bottom frame is for d-wave. Curves are shown 
for T/T c = t = 0.01 (solid curve), 0.3 (dotted), 0.5 (long- 
dashed), 0.7 (dot-short-dashed), 0.9 (dot-long-dashed). The 
short-dashed curve is approximate analytic form for low n 
given in Table II. 
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FIG. 7: The parameter T*/T c versus n for several temper- 
atures shown for the s-wave (top frame) and d-wave (bot- 
tom frame) gaps. Curves are shown for T/T c = t = 0.01 
(solid curve), 0.3 (dotted), 0.5 (long-dashed), 0.7 (dot-short- 
dashed), 0.9 (dot-long-dashed). The short-dashed curve is the 
small n limit (see Table II). T* /T c goes to T/T c as in 0, 
which forms the lower limit on the curves with the upper limit 
being T*/T c = 1 at which point, A(T*) would be zero. 



result at T = derived above. We see that it com- 
pares well with the exact result (solid curve). The other 
curves apply to T/T c = t = 0.3 (dotted), 0.5 (short- 
dashed), 0.7 (long-dashed) and 0.9 (dot-dashed). In this 
case the A(n, T)/A(0, 0) curves do not cross and are all 
constructed from BCS curves for the temperature depen- 
dence of the gap. The temperature T refers to the sample 
temperature before the injection of excess quasiparticles 
n. The intersection of the various curves with the vertical 
axis simply gives the temperature variation of the gap in 
BCS. At finite n, the extra quasiparticles are accommo- 
dated into the system by assuming a higher temperature 
thermal distribution T* , with T* made sufficiently larger 
than T to have n extra thermal quasiparticles. 

Note that in contrast to the /i* model, the differences 
between s- and d-wave are much less pronounced at T = 
0. This reflects the fact that in a thermal distribution, 
blocking effects are not an important consideration. In 
fact now the gap in the d-wave case terminates at a value 
of n which is smaller than in the s-wave case. This is 
opposite to what is found for the [i* model. Also the 
curves show no first order transition to the normal state 
which now occurs only when the gap is zero. 

In Fig. 7 we show the value of T* as a function of 
the noncquilibrium distribution n for various values of 
T. The temperatures used are T/T c = 0.01, 0.3, 0.5, 



0.7, 0.9. Note that the curve with the lowest sample 
temperature (solid curve) at small n agrees well with 
our analytic expressions for the same quantity shown as 
the dashed lines. These fo llow from the transcendental 
equation n = 0.94 v /r*/r c e" 1 - 76Tc/T * for s-wave and the 
explicit equation T* /T c = 2.36n 1 / 2 for d-wave. These 
results are also entered in the final line of Table II. 



III. S-I-N TUNNELING JUNCTION 

Now we consider a specific application of our results 
to the case of a superconducting-insulator-normal metal 
tunneling junction. Denote the current in a S-I-N junc- 
tion with noncquilibrium distribution on the supercon- 
ducting side, described by the /i* model, by I^(V) 
where V is the voltage across the junction. It is given 
by a straightforward modification of the usual tunneling 
formula^O] 

/oo 
deN s (e)[f(e-fj,*)-f(e + V)} (26) 
- OO 
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where iVs(e) is the normalized density of states given by 

* s( ' , = H <^%> (27 » 

with < ... > the average over angles as before. 

We have seen in the previous section that the introduc- 
tion of a nonequilibrium fx* modifies the gap but does not 
change its symmetry and Eqn. (26) still holds for the den- 
sity of states in Eqn. (25) although the new gap ampli- 
tude is reduced by a factor of (1— 2n) and (1— 4\/2« 3 ^ 2 /3) 
for s- and d-wave, respectively, at zero temperature and n 
small. Besides the change in N(e) just described, one of 
the thermal factors in Eqn. (25) is also displaced by the 
new chemical potential fi*. The structure of Eqn. (25) 
makes it useful to separate these two factors, and it is 
convenient to rewrite Ij~*(y) in the form 

/oc 
deN s (e)[f(e-^)-f(e)} 
-00 

/OC 
deN s (e)[f(e)-f(e + V)} (28) 
-oo 

The second term in Eqn. (27) has the identical form that 
applies to an ordinary S-I-N junction in equilibrium at 
temperature T. We denote the current in this case by 
1(V) 

/OC 
deN s (e)[f(e)-f(e+V)} (29) 
-oc 

where the gap amplitude defining N~s(e) is that appro- 
priate to the nonequilibrium superconductor. The first 
term in Eqn. (27) is simply a number, independent of 
voltage. Reference to the defining Eqn. (3) for n shows 
that this number is equal to nA(O). Thus we find 

I™(V) =l(V)+nA(0) (30) 

We see from Eqn. (29) that the current voltage char- 
acteristics are modified in two ways by the nonequilib- 
rium distribution. The entire equivalent equilibrium dis- 
tribution is shifted up by an amount nA(0). This al- 
lows one to measure n once the gap is known. Secondly, 
the "equivalent equilibrium" current voltage characteris- 
tics are those of an equilibrium junction with the smaller 
nonequilibrium gap used instead of its equilibrium value. 
This knowledge allows one to fully characterize the non- 
equilibrium current voltage characteristic and to apply 
checks to see how well the fi* model works. For exam- 
ple, the derivative of I^(V) with V at zero temperature 
simply gives 

dI SN 

-JL-=NS(V), (31) 

the quasiparticle density of states with the nonequilib- 
rium gap but otherwise it is the same as for an equi- 
librium distribution. For s-wave it will have an inverse 



square root singularity at A(n) and for d-wave it will 
go like ln(8A(n)/|A(n) — u)\)/ir instead (see Abanov- 
Chubukov|2lj). In both cases, A(n) can be determined 
from these singularities. Comparison with its equilib- 
rium value gives a measure of n in both s- and d-wave. 
Next it should be possible to check if this is consistent 
with the value of the chemical potential related to n by 
H*/A(0) = l-2n and n*/A(0) = V^n 1 / 2 , respectively, 
for s- and d-wave, at T — and in the limit of small 
n. The chemical potential is measured directly by noting 
that in Eqn. (27) for V = —(J,*, first and second terms on 
the right-hand side are equal but of opposite sign, giv- 
ing a sum of zero. In Fig. 8, we show numerical results 
for IpF(V) versus V at a low temperature T/T c = 0.1. 
The top frame applies to s-wave while the bottom frame 
is for d-wave. It is verified that these curves obey the 
expected rules mentioned above. For the s-wave case 
A(n)/A(0) is set equal to 0.8 while for the d-wave case 
we have used A(n)/A(0) = 0.9 instead. Reference to 
Fig. 1 for s-wave and to Fig. 3 for d-wave shows that 
these choices correspond to an excess quasiparticle num- 
ber of approximately 0.09 and 0.12, respectively. The 
excess quasiparticle number is greater in the d-wave case 
than in s-wave even though the gap is only reduced by 
10% as compared with 20% for s-wave. 

IV. PUMP/PROBE OPTICAL 
MEASUREMENTS OF n(T) 

In the following, we wish to discuss recent experimental 
pump/probe laser experiments which have been used to 
infer information about the excess quasiparticle density. 
In particular, we wish to address a claim that these ex- 
periments provide evidence for s-wave pairing in the high 
T c cuprates. To address this issue, following Kabanov et 
aL0, we use the T* model. While we do not report ex- 
plicitly on this here, we have also examined these proper- 
ties within the //*-model and have found similar results. 

To calculate the excess quasiparticle density n(T), we 
have used Eqn. (3) modified for the T* model via Eqn. (4) 
such that 

n(T) = ^n°°[/(£|,r) - f(E k ,T)]de k \ (32) 



where Et — -J ei + A?(T*) (the asterisk always refer- 
ring to quantities depending on T* instead of T) and an 
average over the angle <f> is done in the case of d-wave. 

To evaluate n(T) at a temperature T, it is necessary 
to know T* and A(T*) and this is determined from the 
amount of laser energy Ei deposited in the system. In 
this work, the laser energy will be assumed to go into 
both electron and phonon systems 

Ei = AE e i ec i ron + AiJp/jonon . (33) 

To begin with, however, we examine the case where the 
energy is assumed to go only into the electronic system: 
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FIG. 8: SIN tunneling I-V characteristics for T = shown for 
the fi* model with s-wave gap symmetry (upper frame) and 
d-wave (lower frame). The current / is normalized by N(0) 
and by the maximum zero-temperature gap in the standard 
way and the voltage V is normalized to the maximum gap at 
T = 0. An excess quasiparticle density n leads to a reduction 
in the gap by A(n)/A(0) and a vertical shift by n in the I- 
V characteristic. I — at V = — . The dotted curve is 
the normal state n = 0, the solid curve is superconducting 
state with n = and the dashed curve is for a reduced gap 
A(n)/A(0) = 0.8 in the s-wave case and 0.9 in the d-wave 
case. 



the quasiparticles and a modification of the supercon- 
ducting condensate due to a A(T*). The energy going 
into the quasiparticles relative to the reference nonequi- 
librium state at temperature T is 



AE„ 



AN{0) 



^J™[Elf{E* h ,T*) - E k f(E k ,T)]de k 



(34) 

Kabanov et treated this piece as [n(T*) — 

n(T)](A(T) + k B T/2) which is not completely correct 
near T c . 

Wc find that the average quasiparticle energy per par- 
ticle calculated as 



E qp = f™ E k f(E k ,T)de k 
N f °° f(E k ,T)de k 



(35) 



gives a constant equal to A(0) at zero temperature for s- 
wave, since excitations can only exist at the gap edge be- 
cause the density of states is zero below this energy. This 
behaviour is seen in Fig. 9 for the dashed curve which 
gives Eqn. (34) normalized to A(0). As T increases, 




FIG. 9: The average quasiparticle energy per particle nor- 
malized to the maximum zero temperature equilibrium gap, 
E qp /NA(0), versus T/T c for s-wave (dashed curve) and d- 
wave (solid curve). The T = value is set by the lowest 
available energy state in the quasiparticle density of states, 
whereas near T c , the energy scale is set by fcsT. 



the energy per particle increases slightly and then de- 
creases near T c to a value of 7r 2 fcsT c /12 ln(2)A(0) = 
1.19(fcsT c /A(0)) ~ 0.67 as now the gap in the quasi- 
particle density of states has shrunk to zero and the en- 
ergy of the quasiparticles is controlled by fc^T which is 
less than A(0). Similar physics is found for a d-wave 
order parameter with the essential difference that exci- 
tations can now occur at zero energy and therefore the 
average quasiparticle energy per particle starts from zero 
at T = and rises linearly reflecting the linear increase 
in energy of the density of states. It can be shown ana- 
lytically that E qp /N A(0) ~ 1.03T/T C for T < A(0), the 
regime where a nodal approximation is valid. At T c , the 
quasiparticle energy per particle is once again controlled 
by ksT and so the limiting number is given by the same 
formula as above but with the BCS d-wave gap ratio of 
A(0)/fc s T c = 2.14 instead of 1.76 for s-wave. The num- 
ber at T c is approximately 0.55 These results are shown 
in Fig. 9 (solid curve) and we will refer back to them at 
a later point. 

The reaction of the condensate is simply given as 



AE r 



2N(0)( 



[Ek - Et 



2Et 



£jjj-(l-2f(E k) T))]de k 



(l-2/(££,T*)) 



(36) 



This term reflects the fact that the presence of excess 
quasiparticles causes a readjustment to the supercon- 
ducting condensate through a change in the gap A* = 
A(T*). This term was not included by Parker [l3j and 
neither was it included in the work of Kabanov et aZ.Q- 
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FIG. 10: The excess quasiparticle fraction, n(T)/no as a func- 
tion of T/T c . These curves are calculated in the T* model for 
fixed laser energy, where the energy is assumed to go only into 
the quasiparticles. The upper frame is for a BCS s-wave gap 
and the bottom frame is for d-wave. Curves are shown for 
fixed laser energy Ei (in units of the condensation energy) 
of 0.2 (solid), 0.1 (dotted), and 0.05 (dashed). Here no is 
n(T = 0) for the case of Ei = 0.2. 



Our procedure was to fix the laser energy Ej and deter- 
mine the T* and A(T*) which gave Ej = AE gp +AE con d- 
For our purposes, we used the BCS temperature depen- 
dence of the gap, calculated numerically, for both A(T*) 
versus T* and A(T) versus T with no approximate form. 
Our results for both s-wave and d-wave gap symmetry 
are shown in Fig. 10 for a variety of Ei, which is nor- 
malized to the zero temperature condensation energy in 
the equilibrium state. Note that the curves shown here 
are normalized to uq = n(T = 0) for the Ej = 0.2 case, 
rather than the n(0) associated with each Ei, in order 
to show the overall relative reduction as Ej is reduced. 
For s-wave, the curves are relatively flat albeit with some 
small depression followed by a sharp upturn near T c and 
then by a drop. The peak occurs when T*/T c = 1, at 
which point we assume that the nonequilibrium state has 
been forced to become a normal metal at an effective 
temperature T* and it is measured relative to the equi- 
librium superconducting state which would exist at tem- 
perature T. Therefore, AE electr on = AE electron (A* = 
0,T*,A(T),T) with T* being fixed by E x . Likewise, 
n{T) = n(A* = 0, T*, A(T), T). The behaviour of the s- 
wave curve largely mimics the inverse of the curve for the 
energy per quasiparticle. At low temperature, the num- 
ber of excess quasiparticles is relatively constant with a 



slight decrease as the temperature is raised, reflecting 
the fact that the energy per quasiparticle is increasing 
slightly and so fewer quasiparticles can be created at 
fixed energy. At high temperature near T c , the energy 
per quasiparticle is decreasing and so more quasiparti- 
cles can be created for fixed energy and one finds that 
n(T)/n(0) shows an upturn in response to this. Likewise, 
the d-wave curve for n(T)/n(0) can be understood from 
the behaviour of the E qp /N curve, with the n(T)/ra(0) 
decaying dramatically as T increases reflecting that it is 
costing on average more energy to create a quasiparticle. 
We can easily show for the d-wave case that 



n(T) 
n(0) 



= 1 



2.9^ 



2/3 



2.9t c 

~e7 



2/3 



(37) 



for small reduced temperature t — T/T c with T <C A(0). 
For t < E It n(T)/n(Q) ~ 1 - (2.9/£/) 2 / 3 i 2 . For t > £j 
but still with T < A(0), n(T)/n(0) ~ (2/M)(E I /2.9) 1 / 3 
(inverse t law). Our numerical results conform to these 
limits. Also note that n(0) = 0.18(£//2.9) 2/3 so that 
n(T) unnormalized to n(0) will go like Ej in the region 
where the 1/t law applies. Once again as the energy scale 
reverts to ksT near T c the slight upturn in n(T)/n(0) is 
reflecting the smaller energy required to create the quasi- 
particles. Kabanov et al. ,7] do not find this result due 
to their approximations and the details of their curve 
would differ as they have only included an approximate 
linear form of the d-wave quasiparticle density of states 
rather than the full form with temperature dependence 
as is done here. In fact, if their data did not go so low 
in temperature and given that n(0) is not known exper- 
imentally, the flatness of the d-wave curve with a slight 
upturn near T c placed on an arbitrary scale, would prob- 
ably make as viable a comparison with their data as the 
s-wave case. However, we note that they do show data 
at lower temperature and so this interpretation does not 
hold, also the 2/3 dependence on Ej at T — is not veri- 
fied. On the face of things, it may appear that their data 
agree best with s-wave. However, we argue, as they did, 
that it is necessary to include phonons in this picture. 

In their analysis, to obtain agreement with their data, 
Kabanov et aZ.pfl did include the fact that some of the 
laser energy would be distributed to phonons in the sys- 
tem. In this case, we partition the laser energy with the 
phonons as well: 



E T = AE n 



AE r 



AE, 



phonon 



(38) 



The phonon piece is calculated assuming that only 
phonons with energy hw above 2A can be considered to 
be out of thermal equilibrium with the lattice and there- 
fore at a temperature T* . It is in this way that a bot- 
tleneck at 2A is introduced into the model, which then 
deviates from pure heating. Such was the same consid- 
eration of Kabanov et al. 7}. As such we calculate the 
amount of energy going into the phonon system as: 



AE, 



phonon 



2 A 



uF(uj)[n{uj,T*) -n(w,T)]duj, (39) 
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where the usual Bose-Einstein factor n{u>, T) — 
1 / (exp(?UJ I ksT) — 1) and F(u>) is the phonon frequency 
distribution function measured by Renker et aL|22l| from 
neutron scattering experiments on YBCO. In our calcu- 
lation, we effectively fix iV(0) to get the correct ratio 
of phononic specific heat at T c relative to the electronic 
part via comparison with the specific heat data of Lo- 
ram et a/.|23l|. This is to ensure that, the phononic and 
electronic portions are balanced in accordance with ex- 
periment. As the phonon energy increases typically as 
T 4 , one sees that this term, as long as T* > T, will take 
more and more of the fixed laser energy away from the 
electronic system and hence, there are fewer excess quasi- 
particles that can be created at higher temperature and 
the curve for n(T)/n(0) must go down. This is illustrated 
in Fig. 11, where curves decay rapidly as T increases and 
are further reduced for lower Ej . Also shown on the same 
figure are the experimental results of Kabanov et ai.0 
(solid squares). We conclude that their data does not 
support an interpretation of s-wave gap symmetry in the 
high T c cuprates. Nor does it agree with d-wave (Fig. 
10, bottom frame). 
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FIG. 11: The excess quasiparticle fraction n(T)/n(0) using 
a T* model in s-wave including phonons for parameters ap- 
propriate to YBCO. For Pb parameters the curves are very 
similar (and not shown) . The curves are for fixed laser energy 
(in units of the condensation energy) of 0.3 (long-dashed), 0.2 
(solid), 0.1 (dotted), and 0.05 (short-dashed). The Pb data 
of Carr et al. |ll|. reproduced here as the black dots, shows 
a suppression with increasing T in keeping with our s-wave 
results with phonons. The YBCO data of Kabanov et 
are shown as the solid squares and disagree with the theory. 

We have also done this calculation with the Pb phonon 
spectrum[24| (adjusted to the specific heat in Pb and 
using a BCS gap ratio) and we find similar curves to 
those shown here. The excess quasiparticle density in 
Pb has been measured by Carr et aL fill and compared 
successfully to rate equation calculations [25| used for de- 
termining the nonequilibrium distribution. We show the 



Pb data (solid circles) on our curves to emphasize that 
Pb, as an s-wave superconductor, does follow the trend 
of showing a suppressed excess quasiparticle density as 
the temperature increases and agrees well with our cal- 
culations. This comparison also serves to show that our 
simple procedure of introducing the bottleneck at 2A and 
sharing the laser energy between phonons and electrons 
agrees qualitatively and even semiquantitatively with the 
more sophisticated and accurate rate equation calcula- 
tions used by Carr et q/.|ll| and validates our simpler 
method. 

Here, we have done the calculation using a BCS 
gap ratio of 3.53. A full strong coupling Eliashberg 
calculation j2|J would have to be done to include a larger 
ratio, as from our experience, simply inserting a larger ra- 
tio in a BCS calculation can give incorrect, and therefore 
misleading, results. Aside from the inherent complex- 
ity of such a calculation, we would need to commit to 
some specific mechanism since phonons are not believed 
to be the source of the high T c . But there is no con- 
sensus on mechanism. To fit experiment, however, Ka- 
banov et aZ.0] phenomenologically increase the value of 
the ratio 2A(0)/fcsT c to about 9. There is, however, no 
rigorous justification for such a procedure and this is our 
main objection to such a fit. To increase the gap ratio, 
it is necessary to increase the ratio of T c /uj\ n in Eliash- 
berg theory 26] , where uj\ n is a particular moment of the 
electron-phonon spectral function which gives the appro- 
priate measure of the average phonon energy involved. 
When this is done, damping effects, entirely left out of 
BCS, become dominant and superconducting properties 
acquire behaviours that are qualitatively different from 
straightforward extrapolations of BCS behavio ur ( see, for 
instance, many properties calculated in Ref. |26j in the 
limit of large T^/wi n ratio). For YBCO, the gap ratio 
is closer to 5|23 and is certainly not 9. Further, for a 
gap ratio ratio of 9-10, the cutoff of 2A applied to the 
phonons falls at 70-80 meV which is at the very top of 
the measured phonon spectrum)^. This large value of 
the cutoff has the effect of greatly reducing the ability 
of the phonons to share in the laser energy and this par- 
tially accounts for why the curve for n(T) in this case 
stays flat to much higher temperature than for the BCS 
curve. 



V. CONCLUSIONS 

In summary, we have examined the differences between 
an s-wave order parameter versus a d-wave in a nonequi- 
librium superconductor using two prominent models in 
the literature, the T* model of Parker and the fj,* 
model of Owen and Scalapino0|- While these models 
may be considered to be somewhat crude, they have the 
virtue of being simple, and accessible in terms of both 
calculation and physical intuition. As a result, one finds 
them still being used by experimental groups to aid in the 
interpretation of their data. With the advent of high T c 
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cuprates and the deeper examination of the issue of order 
parameter symmetry, a re-examination of these models 
allows for the prediction of different power law depen- 
dences on excess quasiparticlc density expected between 
s- and d-wave gaps. Tables I and II summarize such pre- 
dictions. These predictions are grounded in interesting 
physics such as the blocking of states and how the con- 
densate readjusts and as such, they should remain rele- 
vant even within more complicated models. It is hoped 
that the simplicity of our results may inspire further ex- 
perimental and theoretical efforts to examine nonequilib- 
rium phenomena in the presence of an order parameter 
with nodes. 

In addition, we remind the reader of the past discus- 
sions of an inhomogeneous state in the /i*-model and pro- 
vide the prediction that the d-wave state will not be un- 
stable to such a state for the most part. In fact, the 
d-wave state may form a more stable and intuitively in- 
teresting state in which to probe nonequilibrium super- 
conductivity. 

In our work we have specifically addressed two experi- 
ments. SIN tunneling has always been a powerful probe 
of s-wave superconductors and with our work within the 
/x*-model, we show how one may use this experiment to 
measure the model parameters of fi*, n and A(n) in order 
to test the predictions of the theory, both for s-wave and 
d-wave cases. The modern use of STM may provide a 
more attractive avenue for investigating this issue in the 
face of inhomogeneity where the local density of states 
may vary with position within the same material. 

Finally, recent pump/probe experiments in YBCO 



which have been interpreted as providing support for s- 
wave gap symmetry 7] are reconsidered. Within a BCS 
description of the superconducting state and a T* model 
for the nonequilibrium distribution, our calculations, in- 
cluding phonons, do not produce an excess quasiparticlc 
distribution n(T) which is nearly constant in tempera- 
ture with a peak near T c . Rather a quick decay with 
increasing T is found as more of the laser energy is taken 
up by the phonon system. When the explicit case of Pb 
is considered rather than YBCO, the same rapidly de- 
caying characteristic is found and this is in good agree- 
ment with the recent data of Carr et al. [1 1| in this clas- 
sic s-wave superconductor. Our final conclusion is that 
present pump/probe experiments in YBCO cannot be ac- 
counted for by either s- or d-wave gap symmetry and it 
may be necessary to re-examine the interpretation of the 
data in terms of the excess quasiparticle density. In this 
regard, a next step might be to calculate the optical con- 
ductivity itself in the nonequilibrium state so as to make 
a more direct contact with what is measured. 
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